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A method for deriving Friedmann-Robertson- Walker (FRW) solutions developed in Int. J. Mod. 
Phys. D5(1996)71-84, is generalized to account for models with non-minimal coupling between the 
dark energy and the dark matter. New quintessence and phantom (flat) FRW solutions are found. 
Their physical significance is discussed. Additionally, the aforementioned method is modified so 
that, "coincidence free" solutions can be readily derived. Besides, we review some aspects of the 
phantom barrier crossing. In this regard we present a model which is free from the coincidence 
problem and, at the same time, does the crossing of the phantom barrier oj = — 1 at late time. 
QQ ' Finally, we give additional comments on the non predictive properties of scalar field cosmological 

, models with or without energy transfer. 
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Although many alternative models of the universe that take account of the present stage of accelerated expansion 
are being studied, due to their simplicity dark energy (DE) models are perhaps the preferred ones in the bibliography. 



I Notvifithstanding these are plagued by many problems, among which is the so called "cosmic coincidence" (CC) 



problem. Other undesirable features of these models are related with the possibility- already anticipated by the 
observational data- that the DE equation of state (EOS) parameter might be more negative than minus unity (the 
limiting value of the vacuum EOS parameter). This led cosmologists to face the so called "phantom" dark energy, 
i.e., a DE component with negative kinetic energy term. No matter whether or not there are real chances for the dark 
energy (whether vacuum, quintessence or phantom) to exist, the study of such models is very interesting and it could 
hint at more compelling possibilities. 

Perhaps, the simplest realization of the dark energy is a minimally coupled (self-interacting) scalar field. Due 
ps) to their mathematical simplicity, these models- also called as "quintessence" or "decaying vacuum" DE- have been 

I ' intensively and detailed studied over decades. The next step towards a more complicated and, may be, more realistic 

model, is to consider the possibility of an additional non gravitational interaction between the DE and the background 
' cosmic fluid (see reference |1| for a general formalism to describe the coupled quintessence, phantom, non phantom 
, K-essence and Tachyon, etc, and [H for a coherent review). Otherwise the scalar field DE is coupled non- minimally 
ILjl ' to the background fluid (usually cold dark matter) and, consequently, the background particles do not follow the 
geodesies of the metric, i.e, these are coupled to a scalar-tensor metric (see below for details). 

Although experimental tests in the solar system impose severe restrictions on the possibility of non-minimal coupling 
5h I between the DE and ordinary matter fluids (the background) ^ , due to the unknown nature of the dark matter (DM) 
as the major part of the background, it is possible to have additional (non gravitational) interactions between the DE 
component and DM without conflict with the experimental data. Besides, in references [1] and Q it is shown that 
interacting DE models are well consistent with current observational bounds.^ Since there are suggestive arguments 
showing that observational bounds on the "fifth" force do not necessarily imply decoupling of baryons from the 
dark energy Q, then baryons might be considered also as part of the background DM that is interacting with the 
quintessence field. 

Therefore, as it is done in reference 0], we may consider a universal coupling of the quintessence field to all sorts of 
matter (radiation is excluded). Since the arguments given in the appendix of reference Q to explain this possibility 
are also applicable in the cases of interest in the present study, we refer the interested reader to that reference to look 
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^ It should be pointed out, however, that when the stability of DE potentials in quintessence models is considered, the coupling dark 
matter-dark energy is troublesome 0] . Other arguments based on oscillations or exponential blowup of the matter power spectrum that 
is inconsistent with observations, rule out tight coupling Q], 
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for the details. However we want to mention the basic arguments given therein: A possible explanation is through the 
"longitudinal coupling" approach to inhomogeneous perturbations of the model. The longitudinal coupling involves 
energy transfer between matter and quintessence with no momentum transfer to matter, so that no anomalous 
acceleration arises. In consequence, this choice is not affected by observational bounds on "fifth" force exerted on the 
baryons. Other generalizations of the given approach could be considered that do involve an anomalous acceleration of 
the background due to its coupling to quintessence. However, due to the universal nature of the coupling, it could not 
be detected by differential acceleration experiments. Another argument given in f?\ is that, since the coupling chosen 
is of phenomenological nature and its validity is restricted to cosniological scales (it depends on magnitudes that are 
only well defined in that setting), the form of the coupling at smaller scales remains unspecified. The requirements for 
the different couplings that could have a manifestation at these scales are that i) they give the same averaged coupling 
at cosmological scales, and ii) they meet the observational bounds from the local experiments. We complete the afore 
mentioned arguments, by noting that these are applicable even if the coupling is not of phenomenological origin like 
in the present investigation where the kind of coupling chosen is originated in a scalar-tensor theory of gravity. 

Models with non-minimal interaction between the dark energy (whether scalar field modelled or not) and the 
background DM are appealing since the cosmic coincidence problem -why the energy densities of dark matter and 
of dark energy are of the same order precisely at present?- can be avoided or, at least, smoothed out 0, [13, llH ■ In 
the present study we will show that these models are compelling, besides, because they can do the crossing of the 
phantom barrier (w = — 1) with just a single scalar field, a possibility that is incompatible with models with minimal 
interaction between the scalar field and the background fluid. In the latter case the above mentioned crossing is 
possible only if two or more scalar fields are considered [l^, [l^ ■ 

Many exact Fredmann-Robertson- Walker (FRW) solutions have been found in models in which the background 
fluid and the self-interacting scalar field DE have no other interaction than the gravitational one. However, no such 
gallery of solutions exists for the case where the coupling between the background and the scalar field is non-minimal. 
The present paper is partially aimed, precisely, at deriving classes of exact FRW solutions in models where, besides 
the gravitational interaction between the components of the cosmic fluid, additional non-gravitational interaction is 
also considered. To this end we extend a method formerly applied to derive FRW solutions in models with a scalar 
fleld minimally coupled to the background fluid (Ref. J^]) to account for non- minimal coupling also. The classes 
of solutions found comprise both the phantom and the quintessence. For the input and coupling functions chosen, 
self-interaction potentials of the exponential and sinh like form arise. The importance of these kinds of potential is 
noticeable. Exponential potentials (see for instance reference jT6| ) are found in higher-order [l3| or higher-dimensional 
gravity [i^. These arise also in Kaluza-Klein and in string theories and due to non-perturbative effects (gaugino 
condensation [l^, for instance). Their role in cosmology has been investigated, for instance, in [l^, |23|. Sinh- like 
potentials have been studied, mainly, in the context of quintessence models of dark energy and as candidates for the 
dark matter[2lj. To complement the study, in the final part of the paper, we modify the method of reference [l^ so 
that coincidence free solutions could be derived. Another goal of this paper is to study models which not only are 
free of the coincidence problem but, at the same time, do the crossing of the phantom barrier = — 1 at late times. 
We also give additional comments on the non predictive properties of scalar field cosmological models with or without 
energy transfer. 

The rest of the paper has been organized in the following way. In the next section II, we explain the details of 
the model of interacting components of the cosmic fluid we want to explore, including the fleld equations, etc. The 
model includes the possibility to deal with phantom flelds by considering arbitrary sign of the scalar fleld kinetic 
energy. In Section HI the method of reference [l^ (formerly used to generate exact FRW solutions in models without 
interaction between the components of the cosmic fluid), is generalized to account for models with interaction and 
with arbitrary sign of the kinetic energy of the DE. As an illustration the method is applied to the simplest case with 
minimal coupling. Section IV is devoted to deriving of exact solutions in models with non-minimal coupling between 
the DE and the background fluid. In section V we comment in detail on the physical content of solutions found in the 
former section. In section VI, emphasis is made in explaining how to deal with the coincidence problem in the model 
with interaction between the components of the cosmic fluid, subject of the present investigation. For this purpose, 
the method of section HI is modified so that coincidence free solutions can be derived. In section VII we study a 
model that does the crossing of the phantom barrier. Finally, in Section VIII, we summarize the main achievements 
and shortcomings of the paper. 

II. THE MODEL 

We consider the following action that is inspired in a scalar-tensor theory written in the Einstein frame. The matter 
degrees of freedom and the scalar field are coupled through the scalar-tensor metric x{4')~^9ab 
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+x-^WC„,i^l,v^l,x-'9ab)}. (i) 

where e — ±1 (e = —1 for phantom DE, while e = +1 for quintessence), ¥{(/)) is the scalar field self-interaction 
potential, £,„ is the matter Lagrangian is the collective notation for the matter degrees of freedom), and x{4')~'^ 
is the coupling function. 

This action could be considered, instead, as an effective theory, implying additional non-gravitational interaction 
between the components of the cosmic fluid. When the coupling between the scalar field and the matter is minimal 
(no other interaction than the gravitational one) xi't') = Xo = 1- The field equations derivable from the action ([T]) are 

R,.-l9,.R = T^t^+Ti::\ (2) 

where 

Tl£> = eV^0V,0 - \g^u{e{Vc^f + 2^(0)}, (3) 

are the stress-energy tensors for the ordinary matter degrees of freedom and for the self-interacting scalar field (dark 
energy) respectively. These fulfill the following conservation equation: 

V''T(;:)+V'^Ti^)=0, (4) 

so that energy is not separately conserved by each one of the species in the cosmic mixture. Instead, the following 
dynamical equations hold: 

V'Ti^) = Q„ (5) 

where is the interaction term. This is precisely the basic feature of interacting models: there is exchange of energy 
between the components of the cosmic fluid. 
For FRW universe with metric 

ds^^-df + a^l —r+r^dn^], (6) 

where dfl^ = dO^ + sin^Odtp^ , and the spatial curvature k = ±1 or fc = 0, the Friedmann and Raychaudhuri equations 
look like 

3/fc 

3i?^ + ^ - P,n + /90 + A, (7) 

and 



2ij - 2^ = -(p,„ + p„ +P4, + P4,), (8) 

respectively, where — e(j>^/2 + V and Pif, — e(/)^/2 — V. In the last equation we have assumed — (7^ — l)/30, 
where 7^ is the scalar held barotropic parameter, which is related with the equation of state (EOS) parameter uj^ by 
the relationship 7^ = -I- 1. 
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The null component of the conservation equation for the matter degrees of freedom in equation ([5]), can be written 

as 

/9,„ + 3H{pjn + Pm) = Q, (9) 

where the interaction term Q is given by 

Q = Prn^ ^ P„.H[a^^], (10) 

X da 

and the following "reduced" coupling function x = ;Y(37m-4)/2 ^ < 2 - the matter barotropic index) has been 
introduced. Equation ^ with Q given by PH]) can be readily integrated to yield: 

P-m ^ PrnO a'^'''" X, (H) 

where we have considered that the ordinary matter degrees of freedom (the background) are in the form of a barotropic 
perfect fluid, so that pm = (7m — l)pm- 

For the DE one has, instead (null-component of the second equation in ((S])): 

p^ + 3H{p^+P4,) ^~Q, (12) 

or, equivalently; 

(j}[e4> + 3eH(j) + V']= -Q. (13) 

Equations ([7]), (|10p . ((TT|l and ((T3)l represent the basic set of equations of the model of interacting components of 
the cosmic fluid we are about to investigate. In what follows we shall apply a method for deriving new exact solutions 
to this set of equations. 



III. THE METHOD 

In this section we will generalize the method developed in reference [l^ so that we can derive exact FRW solutions 
in models with interaction, such as the model detailed in the former section. To this end, lets assume that the relevant 
functions can be given in terms of the scale factor: = pj^ia), p<f, — P4>{o-), X = x(^)i 6tc. If we introduce the new 
time variable: dij = a~^dt, then equation can be written as 

where we have introduced the input function F = F(a). This function is chosen in such a way that the self-interaction 
potential V can be rewritten as function of the scale factor in the following formfisj: 

Via) = (15) 
Equation (fT4)) can be readily integrated to yield 



a^P(t, 



da[6--p^oa^<-^-'"-^^]+C, (16) 
a da 



where C is an arbitrary integration constant. If we integrate by parts the second term in the right hand side (RHS) 
of (fTB)) . then we are led with the following equation: 



+ p„0= 3f-[2F + 
J a 

(2-7„)p™oa'^2-7„)^]+^, (17) 
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We now introduce the functions 

G(a) = 3iJ^ = p,„ + P0 + A - 3fc/a^ (18) 

and 



L{a) = e{p^-V)^e{p^-^), (19) 

that wiU be useful in what follows. Both functions G{a) and L{a) — 0^/2 are always non-negative: G{a) > 0, 
L(a) > 0. The cosmological constant A can be absorbed into the self-interaction potential V{(p) so we can set A = 
without loss of generality. 

In what follows, for sake of simplicity and unless the contrary is specified, we choose the spatial curvature k — 0, i. 
e., we will explore flat FRW cosmologies. Another magnitude of interest, that will be useful in the future, is the DE 
barotropic parameter: 

After considering equations (jl6p and (|17|). the functions G{a) and L{a) can be written in the following form: 

ay J a 

(2 - 7m)Pmoa'^'-^"^x] + C/a\ (21) 



and 



L{a) - eG{a) ~ ^{F + Pmoa'^'-^-^ x) , (22) 



respectively. 

Exact solutions can be found in the form of quadratures [l| 



Ar^^iVs/^:^, (24) 
J a^a^G{a) 



and 



= 



In equations ([23]) . ([24l) the "±" sign in the RHS means both time directions and, since Einstein's equations (equations 
0, (HI), pn)) and (UH)) are invariant under time inversion then, in what follows, we choose the branch with the 
sign in equation (l23l) (or ((241) '). 

Once the input function F — F(a) and the coupling function x = x(fl) ^^re given, we can find G — G{a) and 
L ~ L{a) through equations (pij) and respectively. Then, by use of equations ([^5]) (or ([M)) ) and (P5)) . one is 
able to find t — t{a) (or = 77(a)) and (p = '/>(a) by direct integration and, by inversion, a = a(t) (or a — 0(77)) and 
(j) = (j)(^t) (or (j) — (f>{ri)) respectively. 

Let us illustrate how the method operates through the study of the simplest example: a model without interaction 
X — Xa — ^ {Q — 0). I. e., the DE and the background fiuid evolve independently. In this case (see equations (pij) 
and 
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and 



G(a) = p^oa-'^" + 4{ / —Fia) + C/6}, (26) 



^(«) = i{6 / ^F{a)^F{a) + C}. (27) 



For simplicity we set C = 0. To chose the input function we assume, additionally^ 



J a s 



(28) 



where s is some arbitrary constant parameter. The equation (I28|) . in particular, is fulfilled if F{a) = B a'^ , where B 
is a constant parameter. Therefore 

G{a) = p,no a-'^- + — a^-^ 
s 

L{a) = e{^^)Ba'-\ (29) 
s 

Taking into account all of the above, the equation (PO)) yields to: 

1, = (30) 

Since, for quintessence < 7^ < 2/3, then 4 < s < 6. For s > 6 (7^ < 0) the DE is a phantom field instead. Let us 
investigate the most general situation when 7^ 7^ (6 — s)/3 = 7,^. In this case the solutions are given by the following 
expressions: 



A ^ .2/3 3t™ 37„i 1 



37m V Pmo 2(37,„ + s - 6) 2 

1 + ,, ,-(^)"^a^--^-% (31) 

2(37m + s - 6) sp„o 



where 2F1 is the hypergeometric function, and 



A(/) = -^^ ^arcsmh[J a 

K y SPmO 

,32, 

so, by inversion of ([32]), one gets a — a{(j)) and, considering equation (jlSp . one finds the form of the potential: 

F(^) = -l/osinh2«[AA0], (33) 

where Vq = B{spmol^BY, q — l/kn and A = ±k/ y^2e/n. Note that the constant parameter q can be written also in 
the form q — —Jdi/ilm ^ 70) ■ Then, if (p were a phantom (7^ < s > 6), g were always positive. 

As explained in this potential is a good quintessencial candidate to be the missing energy in the universe. 
Since it behaves like an inverse power-law potential at early times, then this allows to avoid the fine tuning problem. 
The parameters Vo, A and g, can be determined uniquely by the measured values for the equation of state and of the 
amount of vacuum energy necessary to obtain a tracker solution [2l| . 



^ Other assumptions to fix the input function can also be made. In principle the input function is arbitrary, the only requirement being 
simplicity of the integrals to be taken. 
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IV. GENERATING SOLUTIONS IN NON-MINIMAL COUPLING CASES 



In this section we will apply the method explained in the former section, to generate solutions in models where the 
DE (either quintessence or phantom) and the background fluid share additional non-gravitational interaction through 
the non-minimal coupling given by the coupling function x~'^ (o^'i alternatively, the "reduced" coupling function x) 
in the action ([1]). Through this section, the constant parameters n, fc, etc., are different, in general, from those in the 
former section. 

As before, for the sake of simplicity, we consider C = 0. In this case (see equations (|2T|) and (|22p): 



L(a) = e{G(a) - ^ - p^oa-^^-x}, (34) 



a6 



and 



G(^='^ ^) ^' ^ ^ 

Now we are in position to introduce different input F(a) and coupling functions "xia) to generate solutions. However, 
to illustrate the possibilities offered by the method, it will be enough to choose just one given input function and a 
couple of coupling functions. 

A. F = B a% X = Xo a'^-""" 
After this choice of the input functions F and X: straightforward integration in the first equation in (|34p yields to: 

G[a) = — a'-^ + — ^ a-", (36) 
s 6 — n 

where n 6 and A — 3xo/Omo(2 — 7m). Meanwhile, Eq. psp can be written as follows: 

L _ {6- s)B a"-^ +aa-''^ 



b—n 

s{n - 3jm)A 
rv 

3(6-n)(2-7„)' 
and Eq ([20l) for the scalar field barotropic parameter: 



(37) 



oB ° + a a " 

We point out that, since G{a) > 0, then the constant parameter n is restricted to be n < 6. 



1. General case with s ^ 6 



(i) a = ^ n = 37m 

In this case there is no additional (non-gravitational) interaction between the scalar field and the background fluid: 
X = XOj i-e., this is the simplest situation where minimal coupling between the dark energy and the background fluid 
is considered. According to equation (|38p the DE barotropic parameter is related with the constant parameter s 
through: 37^ = 6 — s. In consequence, for quintessence (0 < 7^ < 2/3) 4 < ,s < 6, meanwhile, for the phantom 
(70 < 0) s > 6. 
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Obviously, if one correctly arranges the constants, this case is the same one considered, as a working example, in 
the last part of the former section. To simplify writing of the solutions, lets introduce the constant k — — 1. The 
integral in (|25p is now easily taken to yield, for k ^ 0: 



or, after inverting it to get a = a{(f>): 



V{(p) = Vo sinh-2/'=[AA0], 

X^±-^^^=^. (40) 
sA 2^e(6 - s) 

Equation (|23p can now be rewritten as follows: 

At = ±-J- :, (41) 

where we have introduced the new variable X — a"/^ and the new constant parameter p = k/(k + 1). The integral 
(|4T|) yields to: 



where, as before, 2^1 is the hypergeometric function, 
(ii) a^O, k = {n = 6- s) 

The "reduced" coupling function x is now of the form x — Xo a'^'''""'"*"^. Since, in this case (see equations ([36)) and 

^. . 6B + A 

G{a) = a °, 

s 

L 6-s + a/B 

the evolution of the scale factor evolution in terms of cosmic time is given by 

a(i) =aoAi2/(6-«), 

a, = [{^-^n^-^r^'-'\ (44) 
meanwhile the self-interaction potential is given by: 

1/(0) = V^) ^ = ±^' (45) 

V Of? 

where Vq = B. Note that there is scaling of the form Pm/ P(t> = const. In the present case, the barotropic parameter 
7^ is given by 



(46) 
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2. Particular case with s — 6 



If in equations P6ll38p one considers s = 6, i.e., V — Vq = B then, one is faced with a situation where the sources of 
gravity are a scalar field without self-interaction potential (only a kinetic energy term present), and the background 
perfect fluid (dark matter, for instance), both "living" in a de Sitter background space-time. It is fixed by the effective 
cosmological constant A = B. The following equations hold: 



G{a) = B + Aa-''/{6-n), 

G '6i?+^a"«' " (6-n)(2-7™)- ^ ' 

As customary, by taking the integral in equation (|23p. one obtains the evolution of the scale factor in cosmic time: 



.(0=a„sinhV"[fVf^*]' ""-[(6^^ 



11/,. 



(48) 



meanwhile, the evolution of the scalar field is given by: 



,2 /ea'(6-n) 



a.csinh[y/^^a-"/^]. (49) 

The scalar field energy density is contributed by the effective cosmological constant"^ and the scalar field kinetic energy 
density: 



P0 = ja-- + B. (50) 

Notice that, the quintessence solution (e — +1) arises whenever a'>0=J>0<n< 3jm (recall that n < 6), meanwhile 
the phantom behavior is displayed once 87™ < n < 6.* However, since a' < for the phantom scalar, there is a time 
tc, such that a^{tc) — aJ^ — —a' /QB ^ p<f,{tc) — 0. For earlier times t < tc, the energy density of the phantom field is 
negative definite. This fact alone, rules out the possibility to describe phantom behavior with the help of the present 
solution, unless the free parameters are chosen in such a way that tc is close enough to the Planck time scale. In this 
case one might argue that the classical theory of gravity is unable to give an appropriated (perhaps semiclassical or 
quantum) description of the cosmic evolution. 

As it will be discussed in the next section, this solution represents an example of a model with transition from 
decelerated into accelerated expansion, where the accelerated phase is driven by the combined effect of the kinetic 
energy density of the scalar field and of the cosmological constant. 

In the next subsection we will investigate cases with the same input function F(a) = B a** and a different coupling 
function x(a). 



B. F = B a% X = Xo a^^" [a"" + a"" 
We write the relevant functions in this case: 

6S' .-6 , A „-n 



G(a) = a 



s 6 ~ n 

B' = B + A/6, A = 3(2-7„)p„oXo, (51) 



^ As already said, this effective cosmological constant can be interpreted, alternatively, as a background vacuum fluid. 
* The possibility n < is ruled out since, in this case, according to l)50|l . the negative energy component of the phantom field increases 
with the expansion. 
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and 



L _ ^ a a'-^ + (3 a 



s b—n 
(6 - S)B' J„rA 



s 



6(2 - 7„) ' 



3(6 - n){2 - 7„i) 

and the barotropic index for the scalar fluid: 

_^ aa^~^ + (3a-" 

(a + i?)a--6+/3a-»- ^^'^> 

We should note that positivity of the energy density (non-negativity of the function G{a)), restricts the parameter 
n to be n < 6. Classes of exact solutions can be easily found for particular relationships among the free parameters. 

1. General case with s ^ 6 

As before, since the input function F{a) = B ^ V ~ B a^^^ . It is also useful to introduce the constant 
/c = n/(6 — s) — 1. For n>6 — s, A:>0 meanwhile, for n < 6 — s, A: < 0. Straightforward integration in (|23p yields to: 

2(s-6)(fc + l) 
k + 1 1 3fc + l Q{n-Q)B' 



Let us try, as before, particular cases where the integral in Eq. ([25|l is easily taken, 
(i) a = 

This choice implies the following relationship involving the constant parameters A, B, s and 7m : 

^3(2-7„0(6-s)^ ^ ^ 

Positivity of the constant B leads to the following restriction on the parameter s: 3(2 — 7„j) < s < 6. After the choice 
of the relationship between the constant parameters made, if one integrates (|25p to find Acf) and then inverts to get 
V = V{(j)), one obtains: 



6(6-n)S^ 1/, {s-erA 
^° ^ ^[ sA ] ' ^ ^ 6./3(6-n) - ^^'^ 

Since n < 6, and since A should be real then, for quintessence (e = +l)/3>0=^ 37 m < n < 6, meanwhile for 
phantom (e = — 1) P < ^ n < 3jm- 

(ii) (3^0^n^3j^ 

The coupling function is given now by 

X = Xo(l + a^-"-+'-^). (57) 



11 



Integration of Eq. ([25)1 yields (recall that Eq. ([54)) . representing part of the solution, is yet valid): 



or, after inversion: 



= ±i(^ — h[^^-^ Y (58) 



y(0) = Vq sinh"2/'=[AA(/)], 

^0 = 5 -^ ^ , A = ±— 59 

^ sA ^ ' 2^eas/B' ^ ' 



We see that, for a phantom field (e = — 1), a < and viceversa. 



2. Particular case with s — Q^V = Vq — B 
We are faced with a scalar field without self-interaction potential. The coupling function looks like: 

X = Xoa'^'" + (60) 

Since pm — Pm,o a~^'^"^ x, the background fluid can be viewed as a mixture of a vacuum fluid with a constant energy 
density plus a perfect fluid with barotropic index n/3. 
The other relevant functions are given as follows: 

G{a) = B' 



6 — n 



6— n 



where now the constant a is always negative since: 

7m ^ 



6(2 -7„0' 



(62) 



while A > and < 7™ < 2. If we integrate ([23)1 and then invert, we obtain the evolution of the scale factor in 
cosmic time: 

o(i) = a„ amli^''['l^At\, 

At the same time, by considering equation (pT|) . straightforward integration of (|25p leads to: 



Acj) = ± A^^!^Z^{in[2a; + b + c + 2^{x + b){x + c)] 



b ^^^^ 2y/cb{x + b){x + c) + cx + b{x + 2c) 
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This last integral is easier to take and, in correspondence, a simpler relationship between 4> and the scale factor a is 
obtained, when one considers the particular case (3 — Q {n — 37m). In this case the coupling function takes the form: 
X = Xo(l + a'^''''"). The scale factor evolves according to Eq. ([63|l as before, however, the integral in ([25|l yields to: 



A4>^±—^ arcsmh[W-^^ — a"'"^]. (65) 

In this case only phantom fields can be accommodated since a is negative. 

Although we have not exhausted all the possibilities offered by the method to generate exact solutions, we think 
the aforementioned solutions illustrate quite well its power. In particular, one can try with other input functions F{a) 
and coupling functions x(o), etc. In the next section we will comment in detail the solutions found in a more physical 
context. 



V. THE PHYSICAL SIGNIFICANCE OF THE MATHEMATICAL SOLUTIONS 



For the subsequent discussion it will be useful to give the expressions that determine several important parameters 
of physical (and observational) importance like the Hubble parameter H or its derivative in respect to the cosmic 
time H (these are useful to determine whether the expansion is accelerated or decelerated) , the dimcnsionless energy 
density parameter for the scalar field fl^ = p^/3H^, etc.^Below we collect the general expressions for them. 

H = ±v/G(a)/3, = 1 - p,„oa~^^" ^, 

ij = G{(7™-7^)r!^-7„,}/2, (66) 
so, in particular, acceleration of the expansion: 



-=H + H^= G'{3(7™ - 7^)^!^ - 37™ + 2}/6, (67) 

is allowed whenever 



S > ^ 3(7™ - j^)n^ > 37„ - 2. (68) 

In this section we will show that it always possible to find a region in (free) parameter space where the given 
solutions are adequate to describe correctly the present cosmological paradigm. We have organized the discussion of 
the solutions previously found, in the same way as these solutions were organized in Section IV. This means that we 
first discuss the solution of subsection A.l (i), then A.l (ii), etc. 

A. F = Ba%x = Xo a'^-""" 
1. General case with s ^ 6 



(i) a = n = 37m 

In this case, as already said, there is no additional interaction between the DE and the DM (x = Xo)- The scalar 
field (DE) barotropic parameter is a constant: 7^ = (6 — s)/3. The expression of the DE dimcnsionless energy density 
is given by: 



Other parameters of physical importance, as for instance, the barotropic parameter of the scalar field fluid 7^, have been defined in 
precedent sections (see equations II38I I. 1 146 I I. II53I I. etc.). 
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= ^ \, r, (69) 

where e = 1 — 3(2 — 'ym)PmoXo/A, is the fraction of DE at the beginning of the expansion (e <C 1), and (3 = 
6(2 — jrn)B/{2 — 70)A. Recall that the solution comprises both quintessence (0 < 7^ < 2/3) and phantom (7^ < 0) 
behavior. Notice that, at late times, the dynamics of the expansion is completely dominated by the scalar field: 
— > 1. Accelerated expansion is allowed once 

or, if one considers that the DM is dust (7™ = 1): Sfi^ > l/(s — 3). Accordingly, this solution allows for transition 
from decelerated into accelerated expansion, as required by the cosmological paradigm that dominates at present. 
Since the DE density is given by: 

= ""'^"^e^ - + ^a^^^'--''^^, (71) 

then, for phantom DE (7^ < => s > 6), the energy density is high at early times and is a decreasing function, until 
a moment when it begins to grow. In particular, at late times, the energy density is a manifestly increasing function 
of the cosmic time, and it becomes infinitely large in a finite amount of proper time, signaling at a Big-Rip type of 
singularity in the future of the cosmic evolution. This is a generic feature of (non-interacting) models of phantom 
energy. 

(ii) Q ^ 0, A; = (n = 6 - s) 

In this case, the coupling function is of the following form: x = xo(i^'^"'~^^~^ , while the dimensionless energy density 
parameter for the DE is a constant. It is given by: 

_ 6B + a 

"^-6STA- ^^^^ 

The DE barotropic parameter is also a constant: 70 = 2((6 — s)B + a)/ {6B + a). In fact, since fl^, = const, flm = 
const, ^ rim/^if, = const, then if the DE is currently dominating the cosmic evolution, it should dominate the entire 
cosmic history, contrary to the evidence from nucleosynthesis about matter dominance in the past. In other words, 
in the present scenario there is no room for transition from matter dominance in the past, to dark energy dominance 
at present, but it is possible to find a bounded region of values for the {A, B, s, a) parameters, where this solution 
describes the accelerated expansion at late time. As easily seen, this solution, with a single exponential potential, is 
not appropriated for the description of the presently accepted cosmological paradigm. This situation is similar when 
we have minimal coupling with a scalar field with single exponential self-interaction potential. 

2. Particular case with s = 6 

This solution is very peculiar since the sources of gravity are a scalar field without self-interaction potential (only a 
scalar field kinetic energy term is present) and the background perfect fiuid, both "living" in a de Sitter background 
space-time. It is fixed by the effective cosmological constant A = B. The DE barotropic parameter is given by: 

meanwhile, the energy density parameter: 

. Ba" + a'/6 

^~ Ba^ + A/{6-n)- ^ ' 
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Worth noting that the DE barotropic parameter is dynamical, i.e., it evolves during the cosmic evolution. It is seen 
that smooth transition from stiff matter in the past, into vacuum fluid in the future, is obtained. The condition to 
get accelerated expansion is translated into the following requirement: 

Notice that, in the present case, transition from decelerated into accelerated expansion is allowed. This means that 
we are faced with a model where the acceleration of the cosmic expansion is driven by the combined action of the 
vacuum (effective cosmological constant) energy density Aq — B, and the kinetic energy density of the scalar field. 
If we neglect the kinetic energy of the scalar field: a' — ^ n = 3jm, then we recover the standard A cold dark 
matter (ACDM) scenario, without interaction between the DM and the background vacuum fluid. This fact provides 
the following interpretation of this particular solution: the interaction between the scalar field and the background 
fluid translates into a change only in the kinetic energy of the scalar field, its potential energy is preserved during the 
exchange. 

As explained in the former section, if the free parameters of this solution are chosen in such a way that tc P4>itc) = 
is close enough to the Planck time scale, the solution could be able to describe phantom behavior also. However, 
unlike a large body of phantom models of DE, the cosmic evolution proceeds without big-rip event in the future. This 
fact is easily explained with the help of equation (|50p for the energy density of the phantom field. Actually, notice 
that the energy density decreases in cosmic time t, contrary to the requirement of an increasing energy density to 
reach a big-rip singularity in a finite amount of cosmic time into the future. 

B. F = B a% X = Xo a^^" [a"" + a""®] 
1. General case with s ^ 6 

(i) a 

This solution is adequate to describe correctly the present cosmological paradigm, through an appropriate choice 
of the free parameters. The barotropic parameter 7^ of the DE, as a function of the scale factor, is given by: 

so that, for n -I- s — 6 > 0, and > 0, the dark energy evolves from being stiff matter in the past, to being vacuum 
energy in the future. The same asymptotic behavior holds for the case when /3 < 0, but the DE behaves like a 
phantom fluid in the meanwhile. 

A very appealing feature of the present solution, is that the cosmic coincidence does not even arise. In fact, choose 
fc>0=>n>6 — s. Lets impose just one more condition on the constant parameters at the beginning of the expansion 
a(0) = 0: 

n^iO) = ^ p^oXo = 7 — , (77) 
— n 

that is compatible with nucleosynthesis constrains on the amount of dark energy at such early times. This restriction 
implies that, in the distant future {a{oo) — 00): 

_ , , {6-n)B-{n + s~6)A 
= (6-n)(B + ^/6) ' 

where, since < < 1: then: 



(79) 
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Therefore, as the cosmic expansion proceeds, the ratio between the background of dark matter and the dark energy 
r = flm/^ij, = Pm/p<p tends to a constant value rg = (1 — ft^{oo)) / fl^{oo) < 1. If the free parameters are correctly 
chosen, it may happen that a state with almost constant ratio r « rg has been already reached and, the future of the 
cosmic evolution is just to asymptotically approach vq. If this is indeed the case, then the "cosmic coincidence" is not 
such a coincidence, but it is the result of the predictable destiny of the cosmic evolution. 

(ii) /? = =^ 71 = 37m 

Here we have a similar situation, with the distinction that the dark energy barotropic parameter is always a constant: 



2a , , 

where a > is for standard DE, meanwhile a < is for the phantom. If one chooses the constant A = 3(2 — 7m)pmoXOj 
then the DE dimcnsionless energy density parameter ll^ evolves from a zero value in the distant past into a value: 



or ^ 6(2 - 7m)i3 - (37m + s - 6)^/3 ^ 
""'^"^^ ^ (2-7™)(6S + A) ^ 

through an appropriate arrange of the free parameters. This means that the coincidence problem can be evaded also 
in this case. This is not surprising, because, as already mentioned, models where there is exchange of energy and/or 
moment between the components of the cosmic mixture are adequate to address this issue. 



2,. Particular case with s = 6=>V^ = Vb = -B 



Finally we have the situation where one is faced with a scalar field without self-interaction potential. As already 
said, in this case only phantom fields can be accommodated since a is negative. 



VI. AVOIDING THE COSMIC COINCIDENCE PROBLEM 



Since models with interaction between the components of the cosmic fiuid are appropriated for handling of the 
cosmic coincidence problem (CCP): why the density parameters of the DM and of the DE are simultaneously non- 
vanishing precisely at present? in order to complement the present analysis, we develop in this section a modification 
of the method previously used to derive FRW solutions, so that we will be able to derive FRW solutions that are free 
of the cosmic coincide. 

Here we follow a procedure similar to that of reference ff\. Let us investigate first, under which conditions the 
aforementioned question does not arise in the situations of interest in the present study. For this purpose it is 
recommended to study the dynamics of the ratio function r: 



r=^ = ^, (82) 
P(t, "0 

in respect to the time variable r = Ina (it is related to the cosmic time through dr = Hdt). The following generic 
evolution equation holds for r: 



r' = ^(^-^) = /(r), (83) 

P<t, Pm P4, 

where the prime denotes derivative in respect to t, and / is an arbitrary function (at least of class C^) of r. One is 
then primarily interested in the equilibrium points of equation ()83|) . i.e., those points Tei at which f{rei) — 0. After 
that, one expands / in the neighborhood of each equilibrium point; r = r^i -\- e^, so that, up to terms linear in the 
perturbations e^: /(r) = {df /dr)r^i€i + 0{ei) e- = {df / dr)r^i€i. This last equation can be integrated to yield the 
evolution of the perturbations in time r: 
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= eoi exp [{df/dr)r,i r] , (84) 
where eoi are arbitrary integration constants. It is seen from (|84[) that, only those perturbations for which: 

{df/dr)r^^ < 0, (85) 

decay with time r, and the corresponding equihbrium point is stable. The necessary condition to evade the CCP is 
then given by the requirement that the point Pm/ Pcp = ^'ei ^ 1 be stable against small linear perturbations of the kind 
explained above. 

If we take into account the conservation equations © and (|12p , and the definition of the interaction term Q given 
in equation pop . then, the function / can be given by the following expression: 



/(r) = r{(lnx)'(r + 1) + 3(70 - >«)}■ (86) 

Note that, for a model without interaction ((In^)' — 0) and with a constant DE barotropic parameter 7^ = 70^0 
(consider, for simphcity, dust-like background fluid so that 7„j = 1); /(r) = 3(7^^0 — 1) and the only (stable) 
equilibrium point is the dark energy dominated solution r = il0 = 1. In consequence the coincidence does arise 
in this case. 

Now we modify the method used before to generate FRW solutions, by considering the necessary condition to avoid 
the CCP. The idea is to give the functions / = /(r) and 7^ = ^(f,{r) as the input functions instead of F and the 
couphng function x- Note, in this regard that, according to 1(551): 



da dr 

dT^ — = — , 87 
a f{r) 



then the scale factor (and any other relevant function) can be written as function of the ratio 



a{r) = aexp J (88) 

where a is an arbitrary integration constant. The next step is to choose a function / such that, at least, one of the 
ruts Te of equation f{r) = 0; obeys that {df /dr)r^ < 0. As explained above, this is the necessary requirement for a 
given solution to be free of the cosmic coincidence. 

Note that, after considering and ([57)l . the equation ([55)1 can be rearranged in the following way: 



d^^X ^ 1 _ 3(70 -7,n) 
dr r{r + l) (r + l)/(r)' ^ ' 



or, after integration: 



where, as before, xo is an arbitrary integration constant, and we have considered that the barotropic index 7^ can be 
given as function of the ratio r: 7^ = l4>(j)- This means that, once f{r) and 70(r) are given as input functions, then, 
by means of equation ([90]) we can find the coupling function x — xi''") (or, after ([88|) . it can be given as function of 
the scale factor: x = x(^))- This means, in turn, that pm (equation (fTT|) ) and = Pm/i" are known functions, i.e., 
the function G(r)) (equation (p^ l is also known. On the other hand, the function L ~ L{r) (equation ([TO)) ) can also 
be written in the following way: 



L{r) = t^^{r)p^{r), (91) 

so, once 70(r) is given as input, L = L(r) is also known. Now, since G and L are both known, then we are able to 
derive exact solutions in quadratures through the integrals (|23l) . and pS)) . that can be rewritten as follows: 
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and 



At = ±V3[ ^^=, (92) 



A0: 



respectively. 

Recall that, once the relevant magnitudes are given as functions of the ratio r, these can be given as functions of 
the scale factor also, through either relationship ([57)1 or 



A. CCP-free Model. 



As it was stated in section IVll it is possible to construct cosmological models which are free from the coincidence 
problem by selecting apropiated functions / = f{r) and 7^ = ^4>{r)- In the case we shall discuss here, we chose: 



/(r) = (3{ro - r), 7^ = 7^,0, (94) 

where (3 and tq are positive reals, meanwhile, the constant 70.0 could be either positive or negative. Recall that 
the solution comprises both quintessence (0 < 70^0 ^ 2/3) and phantom (7^,0 < 0) behavior. Notice that, if one 
substitutes the function /(r) given above, in the evolution equation for the ratio r eq. then the equilibrium 

point = tq is stable. Therefore, the coincidence problem is avoided whenever tq < 1. 

Through ([88|) the relationship between the scale factor and the ratio r can be readily found: 



«W = 7 WP- (95) 

[r - ro)^/P 



Whereas, the coupling function looks like: 



1 

From equation ((TS]), we find the expression for G: 



x{r)=Xo^(^^] , (96) 



G^Go[a{r-ro)r-/''(^^j , (97) 

where Gq — XoPmOj besides, from equation (|9ip the following expression for L is obtained: 



Co7» '°"-:°" ^ . (98) 



2r 

Using Eqs. ([971 [98l [221) and eq. dM]), At and can be computed to give: 



2F1 



A^ = ±^^/^(-o + l)-^/^(.- 



-A/2 



A B -I /i r-Q— r 

2 ' 2 ' 2 ' 1+ro 



A rp— r 



(99) 



where A = 3[7™(ro + 1) + l]/[/3(l + ro)], B = -8(7^,0 - 7m)/[/?(l + ''o)], and 
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FIG. 1: A plot of the energy density parameters Qm (dashed curve) and (thin sohd curve), as well as of the ratio function 
r (thick solid curve) versus redshift, is shown for the model of subsection (I VI A|) . For simplicity background dust {ujm = 0) was 
considered. The values of the free parameters are a = 0.76, (3 — 17 and ro — 0.42. The ratio function r approaches a constant 
value ro at negative redshifts. At present {z = 0) r — 0.43, which means that we are already in the long living matter-scaling 
state. The background energy density dominates the arly stages of the cosmic evolution. At z ~ 0.28 both density parameters 
equate and, since them, the dark energy component dominates the dynamics of the expansion. 



A(t> ^ ±^^1270,0(1 + ro)/3{Vl + ?'o arcsinh(^/f) 



r-oarct^nhUPj-^)}, (100) 
y ro(l + r) 

respectively. The expression of the DE dimcnsionless energy density is given by: 

a+ {ro + IjaP 

The evolution of the energy density parameters fim, ^4>i as well as of the ratio function r versus redshift, is shown 
in the figure[TJ The ratio function r approaches a constant value rp at negative redshifts. At present {z ^ 0) r — 0.43, 
which means that we are already in the long living matter-scaling state. The background energy density dominates 
the early stages of the cosmic evolution. At z ~ 0.28 both density parameters equate and, since them, the dark energy 
component dominates the dynamics of the expansion. At late time, since there is accelerated expansion, the following 
inequality is always satisfied: > (87™ — 2)/(37m — 87^). This solution allows for transition from decelerated into 
accelerated expansion, as required by the well-known cosmological paradigm. 



VII. CROSSING OF THE PHANTOM BARRIER lu ^ -1 



The goal of the present section is to study a model which does the crossing of the phantom barrier at late times. 
To perform the study we use the technique developed in section |Vll 

An incomplete list of previous results on the subject of this section includes the following references: Q (the 
viability requirements on the equation of state and sound speed are discussed) ; realizations of the crossing with an 



extra-dimensional origin, has been discussed in [25|; proposals in the framework of scalar-tensor theories include 
[Hi; [13 single field proposal involving high order derivative operators in the lagrangian); [11] (a model with a 
single dynamical scalar field coupled to an a priori non-dynamical background covariant vector field). In [1^ an 
interacting Chaplygin gas is considered. Standar four-dimensional scalar field models are the most popular options of 
the inventory. 

The impossibility of the occurrence of the transition in traditional single field models [lO] has motivated much 
activity in the construction of two field models that do the job. Examples of explicit constructions can be found in 
[T^ , but perhaps the class of models which have received most attention are quintom cosmologies [H, [ij] . 

In the present section, as already mentioned, we restrict ourselves to the framework of scalar-tensor theories [26|. 
We choose the following input function and DE equation of state (EOS) parameter a;^ = 7^ — 1: 

/W=/3(ro-r), u;^ = u;m + 6'^, (102) 
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FIG. 2: The DE EOS parameter is plotted as a function of the redshift for the the values a = 0.76, (3 — 17, ro — 0.42 and 
S — —0.36. Note that the DE EOS parameter goes from uj^ ~ —0.33 (quintessence) ai z = oo, to uj^ ~ —1.2 (phantom) at 
present. The crossing of the u = — 1 barrier is apparent. 



respectively. In this equation tq, S and /3 are arbitrary constant parameters, and ujm is the state parameter of the 
background fluid. The input function /(r) is the same as in so that the CCP may be avoided in the present 

model as well. 

In figure[2]the behaviour of the DE EOS parameter vs redshift is shown. The free parameters have been appropriately 
chosen (a = 0.76, /3 = 17, ro = 0.42 and 6 = —0.36). Notice that co^ goes from —0.33 (quintessence) at z = oo to —1.2 
(phantom) at negative redshift. It is apparent that the crossing of the uj = —1 barrier actually happens at z ~ 0.2. 
This is always possible if the dark energy and the dark matter have additional non-gravitational interaction (the 
coupling scalar field-background fluid is non- minimal) . The explanation is simple enough: in standard scalar- field DE 
models where the coupling is minimal, the crossing can be achieved by adding new degrees of freedom (extra scalar 
fields). When there is non- minimal coupling instead, there is already an additional degree of freedom that is given by 
the coupling function in the action ([T]). 

Using the same procedure as in the former subsection we find the coupling function: 

X / X 35/(/3ro) 



xM=Xo(^)(— ^) , (103) 



With this coumpling function we are able to derive exact solutions using the integrals (|92p and ([M)) : 



/l2a3\ri-^/2(ro - r)'4/2(,. _ ro)-^/2 

^^-^Vg^^ — .r/^5-2) — ^ 

.Fr(l-f,l + 4,2-f,-^), (104) 



/24((r - 1)(1 + LOr,-,) + (r + 1)(1 + uj,n + 26)) 
(3^1+uj^ + 5) 



^ il+ro)id + r{l+u;m + S)) 

^ 5 + ro{l+.^ + S) n(^'-.'') 

^^^""^ F{y,k)}, (105) 



S + ro(l + (jJm + 6) 



where y = arcsin(-\/l/»' + 1), k — {1 + cum) / {1 + ujm + S) , and / — 1 + ro. In the above equation F{y, k) is the Legendre 
elliptic integral of the first kind, while 11(2:, y, k) is the Legendre elliptic integral of the third kind. 

The dimensionless density parameters have similar behavior as in the former model (figure [1]). In consequence the 
present model is cosnsistent with observational evidence on a matter dominated period in the past and late time DE 
dominance. 

The DM and DE energy densities are given by the following expressions: 
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P™ = Po( — )-^^^'^Xo — T (106) 
r — To \r + l/ \ r J 

P4. = Pm/?' (107) 

At late time (large a), the function rp^ cx oc a"'^'-'''^''/''"^ . This means that, whenever 7 + > 0, the energy 
density is bounded into the future and, consequently, there is not big rip in the future of the cosmic evolution in the 
model under study. 

VIII. CONCLUDING REMARKS 

In this paper we were able to generalize a method for deriving Friedmann-Robertson- Walker (FRW) solutions 
previously developed in [l^, to handle models with non- minimal coupling between the DE and the DM. New classes 
of exact solutions comprising both phantom and quintessence solutions can be readily derived. 

Interacting models of dark energy are useful to account for the coincidence problem. We took advance of this fact 
to extended the method of reference [l^ to derive solutions that are free of the CCP. It has been shown that if an 
adequated input function f(r) and a barotropic parameter are chosen, it is always possible to obtain solutions that 
avoid the coincindense problem. Besides, due to an extra degree of freedom that is related with the coupling function 
in the action ([T]), the crossing of the phantom barrier is indeed possible. Although only two specific models were 
studied in detail, we think it is sufhcient to illustrate how the method developed here allows to construct interacting 
models that avoid the coincidence problem and do the uj = —1 crossing. 

We want to point out that, as shown here, in scalar-tensor models of dark energy where additional non-gravitational 
interaction between the DE and the DM occurs, it is always possible to find a region in parameter space that makes 
a given exact solution to be consistent with observations.^ It is even possible to avoid the CCP and to find solutions 
that are free of the big-rip singularity. For this reason, in agreement with previous claims, we conclude that coupled 
scalar field models are not appropiated to understand the nature of the dark energy at any deeper level [s^l ■ 
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